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A  MODEL  OF  FLEET  DEFENSE  BY  INTERCEPTOR  AIRCRAFT 


DESCRIPTION  OF  THE  ENGAGEMENT 

A  raid  consisting  of  an  arbitrary  number  of  enemy  bombers  and  escorts  is 
assumed  to  proceed  at  a  specified  speed  along  some  straight  path.  Surrounding 
the  bomber  force  is  an  envelope  within  which  penetrating  interceptors  may  attack 
the  bombers  with  missiles.  This  envelope  is  approximated  in  the  two-dimensional 
representation  of  the  model  by  a  circle  whose  radius,  Rg,  depends  upon  the  range 

of  the  interceptor  missiles  and  the  formation  of  the  bombers.  Because  of  the 
greater  range  of  air-to-air  missiles  from  the  head-on  aspect  than  from  tail -on, 
the  center  of  the  envelope  is  displaced  a  distance  d  along  the  raid  direction 
from  the  mean  position  of  the  bombers  (see  figure  1).  Escorts  (one  type  only), 
if  any,  are  positioned  some  distance  (escort  station  distance  E)  from  the  center 
of  this  envelope  so  as  to  make  it  impossible  for  CAP  interceptors  to  attack  the 
bomber  force  without  envountering  return  fire.  It  is  further  assumed  that  the 
defensive  screen  of  the  escorts  is  sufficiently  coordinated  that  an  approaching 
interceptor  becomes  engageable  by  its  proper  share  of  the  escort  force  at  some 
distance  from  the  center  of  the  bomber  envelope.  Similarly,  this  portion  of  the 
escort  force  becomes  engageable  by  the  interceptor  at  another  distance  from  the 
center.  *  The  distances  are  the  sum  of  the  escort  station  distance  and  the 
appropriate  head-on  escort  or  CAP  missile  range,  Rg  or  Rj.  (The  head-on 

range  of  CAP  missiles  against  maneuverable  escorts  will  lilcely  be  less  than  that 
against  bombers. ) 

Following  detection  of  the  raid  by  the  fleet,  CAP  aircra  ft  (of  as  many  as  two 
types)  on  specific  stations  are  assigned  to  the  raid  at  various  times  provided  by 
input.  The  assigned  CAP  are  vectored  along  the  shortest  route  to  the  bomber 
envelope.  Strict  justification  of  this  procedure  requires  that  two  conditions  be 
fully  met:  1)  that  CAP  stations  are  not  so  close  to  the  SAM  zone  boundary  that 
a  straight  path  from  the  station  of  an  engageable  interceptor  to  the  bomber 
envelope  would  pass  within  the  SAM  zone,  and  2)  that  the  point  of  arrival  a .  the 
bomber  envelope  prescribed  by  shortest -route  vectoring  is  outside  the  SAM  zone 
at  the  time  of  arrival.  Astute  CAP  stationing  will  insure  satisfaction  of  the  first 
condition,  and  failure  of  the;  second  condition  has  only  minor  consequences.  If 
the  projected  point  of  arrival  of  the  interceptor  at  the  bomber  envelope  by 
shortest -route  vectoring  lies  within  the  SAM  zone,  the  interceptor  would  ideally 
be  vectored  instead  to  that  point  on  the  bomber  envelope  which  it  would  reach  as  this 
this  point  entered  the  SAM  zone.  The  additional  distance  the  interceptor  would 
travel  to  engagement  in  this  unusual  circumstance  is  at  most  the  radius  of  the 
bomber  envelope,  and  in  most  instances  much  less.  Thus,  engagement  times 
computed  on  the  assumption  of  shortest  -route  vectoring  will  rarely  be  in  error 
by  more  than  one  missile  flight  time. 

*  Partition  of  the  escort  force  into  sections  allocated  to  individually  arriving  inter¬ 
ceptors  is  modeled  by  allowing  each  arriving  interceptor  to  engage,  and  to  be  en¬ 
gaged  by,  the  entire  escort  force.  The  fire  each  aircraft  receives  under  this  some¬ 
what  unrealistic  rule  is  the  same,  on  the  average,  as  would  exist  if  the  opposing 
sides  were  divided  into  engagement  units  with  uniform  force  ratios  which  fight 
independent  battles. 
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FIGURE  1 


After  engagement  has  commenced,  missiles  are  launched  by  each  CAP  and 
enemy  escort  at  vulnerable  targets  until  (a)  the  aircraft  is  destroyed,  (b)  it 
exhausts  its  ordnance,  or  (c)  there  are  no  remaining  engageable  targets.  Indivi¬ 
dual  CAP  aircraft  are  disengaged  when  continued  engagement  with  the  bomber 
force  would  cause  the  CAP  aircraft  to  intrude  into  the  SAM  zone  of  the  fleet. 

The  simulation  stops  when  all  bombers  have  been  destroyed,  or  at  a  specific  time 
which,  although  arbitrary,  must  be  provided  by  input.  This  "stop  time"  can  be 
chosen  to  be  one  of  special  interest,  such  as  that  at  which  the  bomber  force 
reaches  ASM  launch  range.  The  model  allows  CAP  and  escorts  which  have 
exhausted  their  ordnance  to  remain  engaged  in  the  role  of  decoys,  or  to  depart, 
at  the  discretion  of  the  program  user.  By  means  of  target  selection  factors 
specified  by  input,  both  CAP  and  enemy  escorts  can  discriminate  between  target 
types  when  more  than  one  type  is  engageable.  This  feature  is  intended  to  model  the 
judgment  in  target  selection  which  both  sides  can  be  expected  to  employ.  Escorts 
may  fire  preferentially  at  the  more  menacing  type  of  CAP  aircraft  if  it  can  be 
recognized  with  some  degree  of  consistency.  CAP  target  selection  may  show  a 
preference  for  either  escorts  or  bombers,  this  preference  reflecting  in  part  the 
objectives  of  the  CAP  pilot  and  in  part  the  success  of  the  escorts  in  frustrating 
these  objective. 

LIMITATIONS 

Several  facets  of  air-to-air  engagements  are  not  modeled  explicitly.  The 
use  of  deck-launched  interceptors  has  not  been  considered  separately  in  the  model, 
since  such  aircraft  can  readily  be  included  by  computing  their  times  of  arrival 
at  the  periphery  of  the  SAM  zone  via  the  CAP  corridor  from  the  carrier,  and 
treating  this  time  as  an  assignment  time  and  this  point  on  the  periphery  as  a  CAP 
station.  The  effects  of  ECM  can  be  modeled  by  appropriate  adjustments  of  salvo 
rates,  missile  kill  probabilities  and  CAP  aircraft  speed.  Fuel  shortages  are  not 
considered,  it  being  assumed  that  prevailing  CAP  cycling  practice  permits  aircraft 
on  station  to  perform  any  mission  which  might  conceivably  be  required. 

Mixed  ordnance  loads  have  also  not  been  considered.  Where  several  types 
of  ordnance  are  usable  on  a  single  aircraft,  they  may  be  accurately  replaced  by 
a  single  ordnance  type  with  a  salvo  rate  equal  to  the  observed  salvo  rate  for  the 
aircraft  with  a  mixed  load,  and  a  kill  probability  equal  to  the  average  of  the  kill 
probabilities  of  the  individual  ordnance  types  weighted  according  to  the  relative 
frequency  of  use.  A  composite  salvo  of  this  sort  will  give  reliable  simulation 
results  except  in  situations  where  not  all  of  the  loaded  ordnance  types  can  be 
used,  as  will  be  the  case  before  an  approaching  target  reaches  the  maximum 
range  of  the  shortest  range  weapon.  These  periods  should  be  inconsequential 
unless  there  are  large  differences  in  the  maximum  ranges  coupled  with  large 
differences  in  salvo  rate  or  kill  probability. 

Computer  programming  requires  specification  of  the  maximum  sizes  of  the 
CAP  and  escort  forces,  and  the  maximum  number  of  engagement  iterations.  The 
program  for  the  model  as  now  written  sets  ±ese  limits  at  50  CAP,  50  escorts 
and  20  iterations. 


INPUTS 


The  data  required  fay  the  program  is  entered  in  4  groups.  The  first  entry 
group  contains 

fa)  the  initial  random  number. 

Group  I  b)  the  stop  time  (in  minutes),  and 

>~c)  the  desired  number  of  program  iterations. 

This  data  is  supplied  on  a  single  card  in  the  format  (2F10.0,  15).  The  initial 
random  number  is  required  by  the  random  number  computer  routine,  which  uses 
in  the  generation  of  a  random  number  the  value  of  tte  previous  number. 


The  second  entry  group  contains 


Group  II 


d) 

e) 

0 


the  number  of  bombers, 

Rr,  the  range  of  the  center  of  the  bomber  force  at  time  t  =  0, 
(in  miles) 

©r>  the  bearing  of  the  center  of  the  bomber  force  at  time  t  =  0, 
(in  degrees) 

<t>  ,  the  raid  heading  (in  degrees), 

Vr,  raid  speed  (in  knots),  and 

the  distance  from  raid  position  at  t  =  0  to  the  SAM  zone  (in 
miles). 


This  data  is  provided  on  a  single  card  in  the  format  (110,  5F10.0). 


The  next  entry  group  defines  the  combat  capabilities  of  the  escorts  and 
consists  of 

a)  the  number  of  escorts, 

b)  the  number  of  salvos  per  fully -loaded  escort, 

c)  the  salvo  rate  per  escort  (the  reciprocal  of  the  average  time 
in  seconds  between  firings  separated  by  a  damage  assessment 
and,  if  necessary,  acquisition  of  a  new  target), 

r  m  J  d)  the  kill  probability  per  salvo  a  gainst  CAP  type  1, 
roup  e)  the  kill  probability  per  salvo  against  CAP  type  2, 

f)  Rg,  the  escort  missile  range  against  a  target  in  head-on 

aspect,  (in  thousands  of  feet' 

g)  E,  the  escort  station  distant, e  (in  thousands  of  feet),  and 

h)  the  target  selection  factor  '.or  escorts. 


This  data  is  supplied  on  a  single  card  in  the  format  (2110,  6F10.0).  The  second 
item  (b)  is  entered  in  an  array  listing  the  number  of  salvos  remaining  for  each 
escort. 

The  final  entry  group  specifies  the  following  data  for  each  CAP  aircraft: 

a)  the  aircraft  type  (1  or  2), 

b)  t  ,  the  assignment  time  (in  minutes), 

a 

c)  Rj,  the  range  at  assignment  time,  (in  miles) 

d)  0  j,  the  bearing  at  assignment  time  (in  degrees), 

e)  Vp  the  speed  to  engagement  (in  knots), 

f)  Rn,  the  envelope  radius  for  bombers  (in  thousands  of  feet), 

D 

g)  d,  die  envelope  displacement  for  bombers  (in  thousands  of  feet), 

h)  Rp  the  missile  range  against  maneuvering  escorts  in  approxi¬ 
mately  head-on  aspect  (in  thousands  of  feet), 

i)  the  number  of  salvos  fully  loaded, 

j)  the  salvo  rate  (per  second), 

k)  the  kill  probability  per  salvo  against  escorts, 

l)  the  kill  probability  per  salvo  against  bombers,  and 

m)  the  CAP  target  selection  factor. 

This  data  is  supplied  on  a  separate  card  for  each  CAP  aircraft  in  the  format 
(13F6. 0).  A  blank  card  follows  the  CAP  data  deck.  The  data  is  entered  in  the 
array  CAP  (I,  J)  whose  elements  are  defined  in  table  I.  Four  engagement  times, 
CAP  (I,  13),  CAP  (I,  14),  CAP  (I,  15,  and  CAP  (I,  16)  are  calculated  in  the 
program;  status  indicators  CAP  (I,  J),  17<J<22,  are  set  at  the  beginning  of 
each  iteration  of  a  simulation  and  are  modified  as  the  need  arises. 

CALCULATIONS 

The  logical  structure  of  the  program  is  indicated  in  the  flowchart  of 
appendix  A.  Several  calculations  are  performed  in  the  program,  and  those 
whose  nature  is  not  self-evident  are  explained  below.  (Comprehension  of  the 
following  discussion  will  be  facilitated  by  reference  to  figure  1.) 


Calculation  of  CAP  Engagement  Times 


The  position  of  the  center  of  the  bomber  envelope  at  the  time  a  CAP  is 
assigned  to  the  raid  is,  in  a  Cartesian  coordinate  system  (X  axis,  north; 

Y  axis,  west), 


X  =  R„  cos  0_  +  (V  t  +  d)  cos  <$> 
c  r  r  r  a 

Y  =  -  f  R  sin  0  +  (V  t  +  d)  sin  4>  1 
c  I  r  r  '  ra  J 
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(1) 

(2) 


•siere  Rf  is  me  range  and  ©r  the  bearing  of  the  center  of  the  bomber  force  at 

r=o,  *  is  the  raid  beading,  Vr  is  the  raid  speed,  t  is  the  assignment  time,  and 

d  is  the  displacement  of  the  envelope  center  from  the  center  of  the  bomber 
force.  The  coordinates  of  the  assigned  CAP  are 

Xj  =  Rj  cos  0  (3) 


cap  a, 

CAPO, 
cap  a, 
cap  a. 

CAP  (I, 

CAP  (I, 
CAP  (I. 
CAP  (I, 
CAPO, 

cap  a 

CAP  (I, 

cap  a. 
ca?  a 
cap  a, 
cap  a 
capo, 

CAP  (I, 
CAP  (I, 

cap  a, 

CAP  (I, 

cap  a 

CAP  (I, 
CAP  (I, 


TABLE  I 

Array  CAP  (I,  J) 

=  Aircraft  type 
=  Assignment  time 
=  Range  at  assignment  time 
=  Searing  at  assignment  time 
=  Speed  to  Engagement 
=  Envelope  radius  for  bombers 
=  Envelope  displacement 
=  Missile  range  against  escorts 
=  Number  of  salvos  folly  loaded 
))  =  Salvo  rate 

L)  =  Kill  probability  of  salvo  against  escorts 

0  =  Kill  probability  of  salvo  against  bombers 

l)  =  Target  selection  factor  for  CAP 

i)  =  lime  CAP  can  fire  at  escorts 

j)  =  Time  escorts  can  fire  at  CAP 

>)  =  Time  CAP  can  fire  at  bombers 

0  =  lime  CAP  disengages 

l)  =  Operating  (I)/Killed  (0) 

>)  =  Engaged  offensively  with  escorts  (1)/Unengaged  (0) 

))  =  Engaged  defensively  with  escorts  (1)/Unengaged  (0) 

L)  =  Engaged  with  bombers  (1)/Unengaged  (0) 

1)  =  Armed  (1)/Unarmed  (0) 

t)  =  Number  of  salvos  remaining 


Yj  =  -Rj  sin  9j  (4) 

where  Rj  is  the  range  and  8  j  is  the  bearing  of  the  CAP  at  tft.  In  a  frame  of 

reference  which  leaves  the  raid  at  rest,  the  distance  the  CAP  must  travel  to 
the  center  of  the  bomber  envelope  is 

D=^(Xc-Xi)2  +  (Yc-Yi)2  (5) 

The  distance  the  CAP  must  travel  to  engage  the  bombers  is 

°b  =  °  “  V  (6) 

where  Rg  is  the  radius  of  the  envelope  within  which  bomhers  are  vulnerable. 

The  distance  the  CAP  must  travel  to  engage  the  escorts  is 

Qe  =  D-E-Rj,  (7) 

where  E  is  die  escort  station  distance  and  Rj  is  die  range  of  the  CAP  missile 

against  escorts  in  a  head-on  aspect.  The  CAP  becomes  engagable  by  escorts 
after  traveling  a  distance 

Di  =  D-E-RE,  (8) 

where  Rg  is  the  range  of  the  escort  missile.  The  speed  of  the  CAP  in  the 
moving  frame  of  reference,  Vj’,  must  satisfy  the  relation 

(VI’+Vr)2  =  V^,  (9) 

where  Vj  is  the  ground  speed  of  the  CAP.  This  condition  may  be  rewritten 

Vt,2  +  2(Vt’  V  +  VT’  V  )  +  V2  = 

I  lx  rx  Iy  ry  r 

Vj’2  +  2Vj’Vr  [(Xc  -  Xj)  cos  +  -  (Yc-  Yj)  sin  4>  ]/D  +  V2  =  V2  =  V2. 


Solving  for  Vj’,  one  obtains 


Vj*  =  -B  +  y  B2  +  Vj  -  V2  , 
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where: 


B  =  [(Xc  -  Xj)  Vr  cos  <*>  -  (Yc-  Yj)  Vr  sin  *j/D.  (12) 

The  engagement  times  are  thus 

CAP  (I,  14)  =  De/Vj5  (13) 

CAP  (I,  15)  =  Dj/Vj*  (14) 

CAP  (I,  16)  =  D^/Vj*.  (15) 

Calculation  of  CAP  Disengagement  Times 

The  CAP  is  disengaged  when  the  center  of  the  bomber  force  moves  a  distance 
(RB  -  d)  within  the  missiles-free  zone.  This  criterion  for  disengagement  is 

strictly  valid  only  when  the  heading  of  the  raid  is  perpendicular  to  the  missiles-free 

zone  boundary,  a  condition  which  should  usually  be  met,  at  least  approximately. 

Serious  departures  from  normal  entry  can  be  taken  account  of  by  the  use  of  an 

adjusted  boundary.  If  the  distance  from  the  raid  location  at  t  =  O  to  the  SAM 

zone  1 3  D  .the  disengagement  time  is 
s 

CAP  (I,  17)  =  (Dg  +  Rg  -  d)/Vr.  (16) 

Identification  of  Firing  Aircraft 

The  probability  that  a  given  aircraft  fires  a  salvo  is  assumed  to  be  propor¬ 
tional  to  its  salvo  rate,  providing  it  has  an  engagable  target.  This  assumption 
is  the  basis  for  determining  which  side,  and  ultimately  which  individual  aircraft, 
is  responsible  for  a  given  salvo.  The  ratio  of  the  total  salvo  rate  of  ail  engaged 
and  armed  CAP  aircraft  to  the  total  salvo  rate  for  all  engaged  and  armed  air¬ 
craft,  both  CAP  and  escorts,  is  computed.  A  random  fractional  number  is 
obtained,  and  the  salvo  is  attributed  to  CAP  or  escorts  according  to  whether  the 
number  is  less  than  or  greater  than  the  computed  ratio.  If  the  salvo  is  fired 
by  escorts,  this  number  is  used  to  identify  the  firing  aircraft.  If  the  salvo  is 
fired  by  CAP,  a  second  random  number  is  compared  with  the  appropriate  salvo 
rate  ratio  to  determine  the  type  of  CAP  aircraft  responsible  for  the  salvo. 

This  number  is  then  used  to  identify  the  firing  aircraft. 


Where  targets  of  type  1  and  Nj  targets  of  type  2  are  engagable,  die 

ratio  Nj/(Nj  +  crN^)  is  computed,  <j  being  the  target  selection  factor  for  the  firing 

aircraft.  A  random  fractional  number  is  obtained,  and  if  it  is  greater  than  this 
ratio,  a  target  of  type  2  is  selected.  Otherwise  a  target  of  type  1  is  assumed. 
Target  selection  factors  are  evaluated  on  the  supposition  that  CAP  aircraft 
regard  bombers  and  escorts,  and  escorts  regard  CAP  type  1  and  CAP  type  2, 
as  targets  of  types  1  and  2,  respectively. 

Kill  Assessment 

A  random  fractional  number  is  obtained,  and  if  it  does  not  exceed  the 
appropriate  kill  probability  for  the  salvo,  a  kill  is  credited.  Where  the  target 
is  CAP  or  escorts,  the  random  number  is  also  used  to  identify  the  downed 
aircraft. 

Time  Advancement 

Following  the  salvo  whose  effects  are  evaluated  between  points  B  and  K  in 
the  flowchart,  the  time  interval  to  the  next  salvo  is  calculated.  This  delay  is 
computed  on  the  assumption  that  the  probability  that  a  salvo  is  fired  within  any 
small  time  interval  At  is  S  At,  where  S  is  the  total  average  salvo  rate  of  all 
aircraft  which  remain  operating  and  armed  after  the  previous  salvo.  The 
assumption  that  this  probability  depends  upon  time  only  through  the  value  of 
the  salvo  rate  produces  a  Poisson  cumulative  distribution  of  delays  S,  namely 

Pr|  S  <  So  |  =  1  -  e  °.  (17) 

This  assumption  is  almost  certainly  incorrect  for  a  single  aircraft;  following 
the  firing  of  a  salvo,  there  is  a  period  in  which  another  salvo  cannot  be  fired 
by  this  particular  aircraft  either  because  of  the  need  to  assess  damage  or  to 
acquire  another  target.  Thus  for  a  single  aircraft,  extremely  short  delays 
do  not  occur.  Moreover,  extremely  long  delays  should  not  occur,  even  though 
they  are  recognized  as  possible  in  equation  (17).  Equation  (17)  therefore  cannot 
accurately  represent  the  form  of  the  distribution  of  delays  for  a  single  aircraft. 
However,  a  knowledge  of  the  proper  distribution  is  not  needed  to  generate  an 
acceptably  accurate  distribution  of  the  time  required  to  fire  a  given  number  of 
salvos.  Any  distribution  of  delays  for  an  individual  aircraft  which  does  not 
reflect  near  regularity  Ln  the  salvo  interval  will  generate  the  same  distribution 
of  times  required  for  a  given  number  of  salvos  as  the  distribution  of  equation 
(17),  except  at  times  which  are  so  short  as  to  be  comparable  to  the  average 
time  required  to  fire  a  single  salvo.  If  the  number  of  salvos  is  not  small,  errors 
in  the  distribution  at  such  short  “times  will  be  of  little  consequence.  If  the  inter¬ 
val  between  salvos  is  quite  irregular,  as  should  be  the  case  in  a  combat  situation, 


a  distribution  of  required  times  accurate  enough  for  the  purposes  of  the  model 
can  be  expected  after  only  two  salvos.  For  this  reason,  equation  (1)  can  be 
taken  as  the  cumulative  distribution  of  delays  without  risk  of  serious  error. 

The  delay  is  determined  in  the  following  manner.  A  random  fractional  number 
R  is  procured,  and  the  delay  6  is  computed  from  the  relation 

«  =-  [log(l-R)]/S  .  (18) 

The  value  of  6  obtained  in  this  way  has  the  distribution  defined  by  equation 
(17). 

OUTPUT 

After  completion  of  the  desired  number  of  iterations,  the  number  of  escorts, 
bombers,  and  CAP  of  both  types  destroyed,  together  with  the  numbers  of  escorts 
and  CAP  which  have  exhausted  their  ordnance,  is  printed  for  each  iteration. 

The  mean  and  standard  deviation  of  each  of  these  quantities  is  then  computed  and 
printed. 
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PROGRAM  AJRAIR 

DIMENSION  CAP(50 ,23),  SSCARM(50>,  ESC0P<5q>,  nEsRSIPO), 

1  NCDES1  ( 20 ) »  NCDESKSQ  ) «  NCPULl(20?i  NCpUL2C20>.  NESDESI2D ) 

2  NESU\|L<20>,  N: QMQS<20> 

READ  10,  RANONE,  TSTOP*  ITERAT 

10  FORMAT  (2F10,0,I5> 

READ  11,  NB0M8,  RGRAO,  BRGRAD,  HDGRaD,  VRaD,  DlS?SM 

11  FORMAT  <110,  5F1G,Q) 

READ  14,  NESC ,  NbSCSL,  EScSLR*  ESCPKii  EScPk2,  EScmrG,  ESCSTA 
1  ESTSFi 

14  FORMAT  <2110,  6F10,0> 

DO  15  1*1,50 

READ  16,  (CAP <  I , J  ) «  J«l, 13) 

16  FORMAT  <13F6,0) 

V F  ( CAP( l , 5) ,EQ , 0 , )  GO  TO  17 

15  CONTINUE! 

17  NCAPsI-l 
NC  API *0 
NCAP2*0 
MBOMDS*  3 
MESDES*] 

MESUML*3 
MCQES1*3 
MCOES2*  3 
MCPUL1*3 
MCPUL2*3 
SSBQMDS*0 , 

SSESOESsO  , 

SSESJNLsO, 

SSCDES1*0, 

SSCD5S2 sO  i 
SSCP'JLl  *0 , 

SSCP JL2*0 , 

no  15  1 *1 , MCAP 

IF  ( CAP ( 1 , 1 ) t  GT  ,  1 1 )  GO  TO  20 

mcapi*m:a?i+i 

GO  TO  13 

20  NCAP2*nj:AP2  +  1 

18  CONTINUE 

XRAD«RGRAD*C0SF(BHGr AD/57, 29578) 

YRAD«-R3RAD*S1NF(8RGR AD/57, 29578) 

VXRAD*V^AO*COSF<^DGraD/57, 29578) 
VYRA0s-YRAD*SlNF(HCGRAr/57,29578) 

COMPJTE  ENGAGEMENT  AND  DISENGAGEMENT  times 
DO  25  1*1, MCAP 

XC5NAS*i(RAD^<VRAC*CAP(  J  ,2)/60,»CAP<  I,7)/6,076)*COSF(MDGRAD/ 

1  57,29578) 

YCFNASsyRADM vRAr*CAP< l«2)/60,*CAP< I , 7 > /6 , 076 ) *S I NF ( HDGRAD/ 


1  57,2557*) 

xCA?A5a;**t|,3)*eQSflCAPl l, 41/57, 2*57*1 
YCAPA3**CA?(!,3)*S!»iF:CAP{  1,41/57,2*57*1 
Ol5T*S55Tf t (XCE»US«XCA?AS>**l*t¥eSH*S*YeA*ASi**2l 
=i*|(ifCr<*5*XC*P*5)*V«40*CYSS»**«?C*PAS)*V**A|ll/6!S? 
vCiP'*V*.b*SflRir{&*«2*C4PC  lt5)**2*¥»A3**2i 
CISC  !,1S)*CAP(  1,2 >*(DlST-CA*U,4l/4, 07*5/ (VC4*«V/*5t) 

CiOf  »,14»*CAPU,<)*iDlST-(fcScSTA*CA*<  J»i))/$,87S*/(VCAPMV/A0. 1 
Ca3C  l.iSlaCAPU^l^tOlST-JESCSTA^ESSHSSj/i,  §741/4  VCAPHV/60. ) 

?5  riSCI,i7)aiD!STS)'*CCAPt  I*6y*CAP«!e7))/6,0?6)/(VR*8/i{»ll 
CALL  Ri<S£T  CHAMC^E) 
r  r!/)*  C-*k^1  JU  *iCT *uA  A 

^0  1933  ITISiS*!,  iTEHAT 
'•f3^'ta3*<S3«S 

va-wpSc  iT|s«0)*d 

^ta=sciT!'iD5*0 
‘r^Ji'LC  )T!  «5?*0  ' 

*5*C«YsJ 

PSC<3>»*^rSC0L 

TS*E'Cs=5CGL*fcSCSL« 

2*  lsI,*iES£ 

PSCARwc 1 1*1, 

2-  ■»c£S9CS«.C  I  l*3ifcSCSL 
i'iC”£5iC  1 7 1^01*0 
^Cr-£>2  C I T I SJO )  *  0 
^Ccilclt  jTi^0)*P 
'vCs<Js.2t  |T!'«D)sO 
no  2-  I*i,«CAP 
CAP?!, 13>*1« 
rAP{l,l?)aO, 
r*=>c 1,531*0, 

CA  =  { ! , P 1 1 *0  1 
CAP ( C,PP*«1» 

?a  CA»(I,P3)*CAPC!,9) 

71**E*3. 

-:  F _**  C-U-CT  JlHCTifcA  h 

'PTg^'il  H'EXT  ENGAGEMENT  TJHE 
?-  TP''P*I,=9 

n0  31  I  si , VC AP 

IP  { VES^C 1 rU ,  0  )  GO  TO  32 

IF  (CAP(1, 14)»CAWU,18)*CAPU, 221  ,i,£. TIME)  GO  TO  35 
IF  (C*P(I, 14), GE?TEMP)  GO  TO  35 
TE-PsC*3{  1 ,2,4) 

j&  IF  ('E3;0-tt9,0)  oC  TO  32 

IF  <9A={I,13>*CAPC!,18>,LE.T!H6>  GO  TO  32 
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1  v  **r.t 


IF  <CAP(I,l5),G£*rEMP)  GO  TO  32 
T£MP*CA3.(  1,15) 

3?  IF  ( CAP ( ! , 16) *CAp ( I , 18)*CAP(  l , 22 ) , L£ , T 1  ME )  GO  To  3‘J 
IF  (CAP( 1,16) ,GE;TeMP)  GO  TO  30  ' 

TE’-'PsC  A  5  (  1.16) 

3 !•  COMT  I NU r 
TI*E*TE<lP 

IF  (TIm:,3T,TSTOP)  GO  To  100 
HO  40  Isl.NCAP 
no  40  Jsl4 , 16 

IF  (CAP(  I,J)  IKE)  60  TO  40 

CAP( I. J*5)sl, 

40  CO'iTIMU: 

r  FLOW  CHART  JUNCTION)  C 

r  ngTERMIviE  WhlCH  SIDE  SCOOTS 

41  MTEST*G 
CPEOLl*  3 . 

CPPOt2sD  . 

IF  (sES:O,GT,0)  GO  TO  46 

00  45  !sl»feCAP 

IF  (CAP{ 1,1) ,GT,1)  GO  TO  43 

CP5OL1*3P£0L1  +  CAP(!,21)*CAP(!#28)*0aPU»18> 

shca31s:a3u,io)*cpeoli 

GO  TO  45 

43  CPP0L.2  =  :Pc0L2*CAP(I,2l)*CAP(  I#22>*CAP(l,l8) 

$R0AS2a:A3(  If10)*CPEni_2 
4C  COvT  I N'J- 
GO  TO  53 

4*  DO  4°  Isl.NCAP 

IF  (2A»(i*20)*CAP(I,18),EO,t,)  «TEST*1 
IF  ( C A* { t  # 1 ) , GT  » 1 . )  GO  TO  47 

CPE0Lls:P50Ul^{CAPU.l9)*CAP(I,21)-CAP(I,l9)*CAP(l,21)  )*CaP(  1,22)* 
1  CAPU.16) 

SRCAP1s:a»( I,10)*CPE0L1 
GO  TO  45 

4  7  CPE0>2s:P“0L2=>(CAP(  1 ,19)*CAP(  i,21)-CAP(  I  »i9)*CAP(  1 ,21 )  )  *CAP ( I »  22 ) * 
1  C  AP ( 1 , 18 ) 

SRCA=,'2s:A3(I<10)*CPEOL2 
4 A  CO^lNUr 

IF  (Vi6S:Oto,fcG.O)  GO  TO  50 
IF  (  ViTSsT.EQ.Q  )  GC  TO  50 
TSRsSRCAPl+SRCAP^+TSRESC 
IF  ( mOE.PT, 60.1)  GO  TO  97 
TESTCEs(SRCAPl+SPCAP2)/TSR 
CALL  WA'jlMG(X) 
tr  (X.GT.T6STCE)  75,52 
r  FlOw  CwaRT  JUNCTION  0 
5.J  TS«sSRCAPl  +  SRCAP<* 
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51  IF  (VIDE. AY, 60,1)  GO  TO  97 
CAP  SHOOTS 

52  TESTl2*SRCAPl/(SRCAPi*SRCAP2t 
CALL  RAVLM3(X) 

IF  ( X , GT i TEST12 )  GO  TO  56 
CAP  TYPsi  1  SHOOTS 
NZ0N6«( :PE0L1*X/T8ST12)*1, 

KsO 

00  55  Mill NCAP 

IF  (5AP(M#l).NE,t.)  GO  TO  55 
IF  { viESCO'EQ.  9  )  GO  TO  53 

J=(CAP(  H'il9)*CAP{Hi21)-C»-P{M«19)*CAP(M*2l))*CAP(H,18>*CAP(H,22) 
GO  TO  54 

53  jrCA9(H,2n*CAP(Kil8)*CAP(H*22> 

54  KsK  +  J 

IF  ( K  , E3 ,  NZONE )  GO  TO  60 

55  COV-'TtNUi! 

CAP-  TYPsi  2  SHOOTS 

56  \'ZONE*(  :P50L2*(1!-X)/(1,*TEST12>  >  +  l. 

KsO 

00  57  Htlj  NCAP 

IF  { CAP ( H, 1 ) , NE , 2 . )  GO  TO  57 
IF  (MESCO.EQ.O)  go  TO  58 

J=(CAP( 4«i9}±CAP(M,21)-CAp<M,19>*CAp(M*21)>*CAP<Nsi8>*CAP(H,225 

GO  TO  59 

5R  JsCAP<M,21)*CAP(M,18>*CaP(M*22> 

59  KrK* J 

IF  (<,E O',  NZONE)  GO  TO  60 

57  COMT  I N  J  El 

6)  CAPtH.PJl^CAPtMiZS)*!, 

IF  ( CAP ( H, 23 ) , NE , Q  « )  GO  TO  63 

CAP { H , ) iO , 

IF  <CAP{M,l),GT,i,)  GO  TO  61 
MGPULlI  lT;ND}swCPtJLl<  IT  I ND > *1 
GO  TO  6? 

5 1.  NCPUu2(  lTIND)sNCPUL2UTlND>*i 
6?  CO'MTINlJ-i 

0 1 SEvjGA 36  EMPTY  CAP  WITH  NEXT  CARD  (500  ) 

500  CA°<H,19')«0  t 

F|_nw  CHART  JUNCTION  E 
CAP  SHOOTS  AT  BOMBER  OR  ESCORT 
6  ^  IF  ( CAP ( M *  19 ) , EQ i 0  * )  GO  TO  70 
IF  ( VESOO.EQ, 0  )  0  TO  70 

IF  (CAP{M,21) , EU •  0 • )  GO  TO  65 
ROM80»NS0MBO 
ESCQ«NE$CO 

TESTE8s30HB0/(B0MB0+CAP(M,13)*ESC0) 

CALL  RANUMB(X) 

IF  ( X , Lrb  TEST6B  )  GO  TO  70 


C  CAP  shoots  at  escort 

65  CALL  RA\jUMB(X) 

IF  (X.5T.:aP(M,1D)  GO  TO  90 
M2ON=s(iSC0*X/CAP(H1ll) )*1, 

K  =  -J 

no  66  Nsl  *  MfcSC 

JsFSCOB(V) 

KrK  + J 

IF  (<,F3.  '(ZONE)  GO  TO  67 

66  CONTINUE 

67  ESCO3(N)s0, 

\!ESC"sN=SCO»l 

meBCI'LsMSSCOL-I 

MEADES ( lTINO>«h£SDES(lTlND)*l 
TSPE5C=TSRFSC-ESCSLR«ESCARH(N) 

GO  TO  R3 

C  FLOW  CHART  JUNCTION  F 

r  CA3  shoots  at  bobber 

7P  CALL  RA>jUM8(X) 

IF  ( X , G T , CAP ( M 1 12 ) )  GO  TO  90 

NBOMRO=\|B3P80-1 

NBOMPSC  lTiND)=weCKOS(  niND)**-! 

IF'  (v80«lBa,EQ.O)  1 U  0  *  90 
C  FLOW  CHART  JUNCTION  G 

C  ESCORT  5H30TS 

75  FSCO'.s^tSCOL 
CALL  RA^IHB(X) 

NZOwEsC  =SCOL*X)+l, 
k  =  U 

CO  76  L=1.NESC 
J=FSOAR«1(w)*ESCOF(l  ) 

KsK  +  J 

IF  (  <  ,  S  3  » NZONfc  >  GO  TO  77 

76  CONTINUE- 

77  NE5RSLI  .i)=NfeSRSL(L)-l 

IF  {VESRSl(L),N£;c)  GO  TO  78 
ESCARMI  J>*0* 

N’ESUVL(  I  T I  MO )  sNES'JM.  {  IT  IND)*1 

TSPE?C=r5^bSC-ESCSLR 

Mg«?C'’l='iE5COL-l 

c  niSE'iGASE  UNARMED  ESCORT  with  FOLLOWING  Cards  (6U0»6'1) 

*00  ESOO° ( L  )  5 0 i 
601  NESCO=M=SCO-1 
C  FLOW  CMfRT  JUNCTION  H 

C  ESCORT  SHOOTS  AT  CAP  TYpE  i  OR  2 

7 8  CApDEls  3  » 

CAPDE2S3. 

DO  79  I*1,NCAH 

IF  ( CAP ( I # 1 ) *  GT  1 1 • )  GO  TO  80 
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CAf»DFl«;*i»n£l*CAFi  I,20UCAP(  1,18) 

no  t r*  7? 

61-  rAPD?223APGfe2*CAP(  1,20  )*CA^(  1,18 i 
7°  COk'T 1  NJ£i 

TESTj2»:AaD6l/(CAPDEl*bSTSF*CAPUE2) 

CALL  R  A  \JUMB  ( X ) 

IF  ( X,QT. TEST12)  GO  TO  85 
FSCORT  SHOOTS  AT  CAP  TYpE  1 
CALL  RAvl'MBCX) 

IF  ( Y , a T i SSCPK1 }  GC  TO  90 
nZOn5«(  :APD£l*X/ESCPKi  Wl, 

**U 

00  82  Ial,NCAP 

IF  ( CAP (  1 , 1 }  , NE  ,  1 ,  )  GO  TO  82 
JsCAP(I,20)*CAP(I,18) 

K  =  K*J 

IF  ( < , 6  J « YZONfc  )  GO  TO  83 
COVTIN'Ji- 
3S  CAn( ! , 1 8 ' *0  « 

’!C”ES1C  iTI\>d)*NCOcSK  ITiNq)*! 

GO  TO  93 

ESCORT  SHOOTS  AT  CAP  TYpE  2 
At  CALL  RA\|IH8(X) 

IF  Of,GT.sSCP*2)  GC  To  90 
M70N?3<:A30E2*X/ESCPK2)*1, 

KsO 

no  86  Itl,NCAP 

IF  (0A°(1 ,1)  ,  'iEf2«  )  GO  TO  86 
JsCA3(I ,20>*CAP< 1,18) 

Ks-f*  J 

IF  (<tE3.  ^ON'E)  GO  TO  87 
*j4  CON,T  I  'Y  J  3 

8  7  C AP ( ! , 1 3 ) *0  , 

Mcn6S2(  lTINn)sNCOE32UTlND)*l 
flow  chart  junction  k 
targets  available  to  either  side 

9i.  0-0  9i  IslfMCAP 

IF  (0AR(I,2l)*-CAP(I,22)*CAPn,18).\'E,0.)  GO  TO  94 
IF  (''ESOO.EQ.Q  >  GO  TO  91 

IF  (0ARn#19)*CAH(I,22>*CAP(!,l8),ME«0t)  GO  TO  94 
IF  (VES:ute,EQ,0)  GO  To  91 
IF  (CAP(I,20>*CAP(I,13>.NE.0.>  GO  To  94 
91  C  0  M  TJ  N 1 !  r 
75  00  TO  29 

ADYamCE  CwOCK,  ADD  CAP  ENTRIES 

9  4  NOEL  A Ys  l 

on  TO  4 J, 

FLOW  CHART  JUNCTION  L 
97  Malays} 
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n  A I  L  RAijUMB(X) 
nfcLAVs-(L3GFIl-X) )/TSR 
TI^ErTHE  +  OELAY 
IF  (  T  I  M5i,  3T  ,  TSTGP  )  GO  TO  100 
HO  98  I  si >  NC AP 

nn  93  J si 4  1 16 

IF  (CAPfl.JJ.LE, (TIME-DELAY))  GO  TO  98 
IF  (OA3( I , J)  ,GT,TIME)  DO  TO  98 
CAP( !, .1*5)31. 

98  COwTINU=. 

HO  99  Isl,NCAP 

IF  (CAD( I ,17) ,GT;i Irtfc)  GO  TO  99 
CAP(I,13)sO. 

99  CO'TIMUr 
GO  TO  <1 

C  FL^w  C^aRT  JUNCTION  M 

IOC  MBnMi'Ss*<0OMi)‘j^NdCMDS(  IT  IND) 

MES0  =  Ss 'iFSOES+NE^DES <  IT  IND) 
MtGUVLs>lE3UNL-*-Nfc3uNL(  IT  IND) 
MCnESl*«CDESl*NCDfcSlt ITJND) 

MCnE32s'iCDES2  +  NCkES2(  IT  IND) 
MCPULls'iCPlJU+NCEuLlC  1TJND) 
hCPUl2s-;CPUL2+NCPUL2UTIND) 

YsNBOMOS  (  I  T  I ND ) **  2 
SSP0mDSsS5H0MUS+‘ 

Y  =  \'EFDE5(  ITIND)**2 
SSFSnE!?sS3ESDES*Y 
Ys-iEFUM,i(  I T I  NO ) *  *2 
SSESJNLsSSESUnL+Y 
Ys-'COESU  ITIND)**2 

SSCI)  =  SlsSSCDE81  +  T 

YsYCnt??( ITIND)**2 
SSCDES?sS3CUES?*Y 
YsNC  =  ULl( I T l NO  )  **2 
SSCP  -LlsSSCPbLl+T 
Ys'.'C3UL  2  (  l  T  I  NO )  **2 
SSCPJL?sSSCPUL2*Y 
1000  CONTINUE- 

C  CO^P-'TE  AVERAGES  AND  STANDARD  DEVIATIONS 

ABnMDSs^aO'4OS/ITEKAT 
AESOFSs'lESDfcS/ITERAT 
AFSU''Ls*65UNL/lTERAT 
ACDESIsnCDESI/ iterat 
ACr>ES2s'iC0FS2/ITFRAT 
AC®ULl*-<C3Ull/lTt£RAT 
ACDUL2=^CPUL?/ITERAT 
XsITERAT 

SDRO^OSsSORTF  (  (  SS60m!)S»X*A8QMDS  )  /  {  X-l. ,  )  ) 
SDfSDES«S3RTF(  {SibSDES“X*AESDF3)/(!<oi,  )  } 
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SDESDNL*S3RTF(  (SSESUNL-X*AESUNU/(*-li  >> 

SDCQESi*S3RTF((SSCDESl*X*ACDESi)/{X.l,S ) 

SDr.D?S2«S3RTF(  {SSCDES2-X*ACDES2)/{X«1,  > ) 

S0nPULl*S3RTF( (SSCPgLl-X*ACPULi}/(X.i,>) 

SO'*.?  JL2  *S3RTF  {  ( S»CPUL2-X*ACPUL2  )  /  ( X»1 , >  > 

PRINT  1351 

1051  FORMAT  (69H  BOMQERS  CAP  TYPE  i  CAP  TyPE  2 

1  ESCORTS ) 

PRINT  1352,  N80M6 ,  NCAPl,  NCAP2,  NESC 

1052  FORMAT  (7X,I5,9x;i6,14X,I6.15X,l5) 

°R I  NT  1353 

in53  FORMAT  (75HQRUN  DESTROYED  DESTROYED  UNLOADED  DESTROYED  UNLOADED 
1  RESTATED  UNLOADED) 

HO  1155  !*1«!TERM 

'sr  I  nT  1354,  I,  N80HDS(I>,  NCDESIU),  NCPUlI(I),  NC0ES2  ( I ) « 

1  nCPU.2  ( I )  *  NESDESU),  NESUNL(.I) 

4  FORMAT  (tX,l2»3X*l5,6X.I5,4X,I5i6X,l5,4X*l5,6X,l3,4X»I5) 

5  CO^TINUS 
PRINT  1356 

1055  FORMAT  (40H0  AVERAGES) 

PRINT  1357,  AdOMBS,  ACOESl,  ACPULl*  ACDES2,  ACPUL2,  A6SDES,  AESU^l 
105  7  FORMAT  (SX,E7,l,1!X,F7ll,2X,F7,l»4X»F7ll,2X,F7tl,4x,F7,l,2X,F7ll) 
PRINT  1359 

105*1  FORMAT  ( 46H0  STANDARD  OBVIATIONS) 

PRINT  1357,  SD8DMDS,  SDCDESl,  SDCRUlI,  SDcDES2,  SdCPUL2,  SDESDES, 

i  soesjN* 

IrNT 
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